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COMPRESSIBLE MHD TURBULENCE: MODE COUPLING, 
ANISOTROPIES AND SCALINGS 
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ABSTRACT 

Compressible turbulence, especially the magnetized version of it, traditionally has a bad reputation with re- 
searchers. However, recent progress in theoretical understanding of incompressible MHD as well as that in 
computational capabilities enabled researchers to obtain scaling relations for compressible MHD turbulence. 
We discuss scalings of Alfven, fast, and slow modes in both magnetically dominated (low 0) and gas pressure 
dominated (high (3) plasmas. We also show that the new regime of MHD turbulence below viscous cutoff 
reported earlier for incompressible flows persists for compressible turbulence. Our recent results show that 
this leads to density fluctuations. New understanding of MHD turbulence is likely to influence many key 
astrophysical problems. 
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1. ISM AND MHD TURBULENCE 

The interstellar medium (ISM) in spiral galaxies 
is crucial for determining the galactic evolution and 
the history of star formation. The ISM is turbulent 
on scales ranging from AUs to kpc (Armstrong et al 
1995; Stanimirovic & Lazarian 2001; Deshpande et 
al 2000), with an embedded magnetic field that in- 
fluences almost all of its properties. This turbulence 
holds the key to many astrophysical processes (star 
formation, the phases, spatial distribution, heating 
of the ISM, and others). 

Present codes can produce simulations that re- 
semble observations (e.g., Vazquez- Semadeni et al 
2000, Vazquez-Semadeni 2000). To what extent do 
these results reflect reality? A meaningful numer- 
ical representation of the ISM requires some basic 
non-dimensional combinations of the physical pa- 
rameters of the simulation to be similar to those of 
the real ISM. One such parameter is the "Reynolds 
number" , Re, the ratio of the eddy turnover time 
of a parcel of gas to the time required for viscous 
forces to slow it appreciably. A similar parameter, 
the "magnetic Reynolds number", Rm, is the ratio 
of the eddy turnover time to magnetic field decay 
time. The properties of the flows on all scales depend 
on Re and Rm. Flows with Re < 100 are laminar; 
chaotic structures develop gradually as Re increases, 
and those with Re ~ 10 3 are appreciably less chaotic 
than those with Re ~ 10 7 . Observed features such as 



star forming clouds are very chaotic with Re > 10 8 
and Rm > 10 16 . The currently available 3D simu- 
lations for a 512 3 grid can have Re and Rm up to 
~ 6000 and are limited by their grid sizes. It should 
be kept in mind that while low-resolution observa- 
tions show true large-scale features, low-resolution 
numerics may produce a completely incorrect phys- 
ical picture. 

How feasible is it, then, to strive to understand 
the complex microphysics of astrophysical MHD tur- 
bulence? Substantial progress in this direction is 
possible by means of "scaling laws" , or analytical 
relations between non-dimensional combinations of 
physical quantities that allow a prediction of the mo- 
tions over a wide range of Re. Even with its limited 
resolution, numerical simulation is a great tool to 
test scaling laws. 

In spite of its complexity, turbulent cascade is 
remarkably self-similar. The physical variables are 
proportional to simple powers of the eddy size over 
a large range of sizes, leading to scaling laws ex- 
pressing the dependence of certain non-dimensional 
combinations of physical variables on the eddy size. 
Robust scaling relations can predict turbulent prop- 
erties on the whole range of scales, including those 
that no large-scale numerical simulation can hope to 
resolve. These scaling laws are extremely important 
for obtaining insights of processes on the small scales. 

By using scaling arguments, Goldreich & Srid- 
har (1995) made ingenious predictions regarding rel- 
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ative motions parallel and perpendicular to magnetic 
field B for Alfvenic turbulence. These relations have 
been confirmed numerically (Cho & Vishniac 2000; 
Maron & Goldreich 200f ; Cho, Lazarian & Vishniac 
2002a; see also Cho, Lazarian & Vishniac 2002b, 
henceforth CLV02 for a review); they are in good 
agreement with observed and inferred astrophysical 
spectra (CLV02). A remarkable fact revealed in Cho, 
Lazarian & Vishniac (2002a) is that fluid motions 
perpendicular to B are identical to hydrodynamic 
motions. This provides an essential physical insight 
and explains why in some respects MHD turbulence 
and hydrodynamic turbulence are similar, while in 
other respects they are different. 

The Goldreich & Sridhar (1995, henceforth 
GS95) model considered incompressible MHD, but 
the real ISM is highly compressible. This motivated 
further studies of the compressible mode scalings. 

2. EFFECTS OF COMPRESSIBILITY 

Three types of waves exist in compressible mag- 
netized plasma. They are AlfVen, slow and fast 
waves. Turbulence is a highly non-linear phe- 
nomenon and it has been thought that it may not 
be productive to talk about different types of per- 
turbations or modes in compressible media as every- 
thing is messed up by strong interactions. Although 
the statements of this sort dominate the MHD liter- 
ature one may question how valid they are. If, for 
instance, one particular type of perturbations cas- 
cades to small scales faster than it interacts with the 
other types, its cascade should proceed on its own. 
Within the GS95 model the Alfvenic perturbations 
cascade to small scales over their period, while the 
other non-linear interactions are longer. Therefore 
we would expect that the GS95 scaling for Alfven 
modes should be present in compressible turbulence 
as well. 

Some hints about effects of compressibility can 
be inferred from Goldreich & Sridhar (1995) seminal 
paper. A more focused theoretical analysis was done 
in Lithwick & Goldreich (2001) paper which deals 
with properties of pressure dominated plasma, i.e. 
in high j3 regime, where /3 = P gas /Pmag > 1- In- 
compressible regime formally corresponds to [3 — > oo 
and therefore it is natural to expect that for [3 3> 1 
the turbulence picture proposed in GS95 would per- 
sist. The only actual difference between subsonic 
turbulence in high (3 and incompressible regimes is 
that for finite (3 a new type of motions, fast modes, 
is present. Those modes for pressure dominated en- 
vironments are analogous to sound waves which are 
marginally affected by magnetic field. Therefore it 



is natural to expect that they will not distort the 
GS95 picture of MHD cascade. Lithwick & Goldre- 
ich (2001) also speculated that for low (3 plasmas the 
GS95 picture may still be applicable. 

A systematic study of the compressibility in low 
(3 plasmas was done by Cho & Lazarian (2002, 
henceforth CL02). First of all, we tested the cou- 
pling of compressible and incompressible modes. If 
Alfvenic modes produce a copious amount of com- 
pressible modes, the whole picture of independent 
Alfvenic turbulence fails. However, the calculations 
(see Fig [l]) show that the amount of energy drained 
into compressible motions is negligible, provided that 
the external magnetic field is sufficiently strong. 

When the decoupling of the modes was proved, it 
became meaningful to talk about separate cascades. 
A remarkable feature of the GS95 picture is that tur- 
bulent cascade of Alfven waves happens over just 
one wave period. Therefore the non-linear interac- 
tions with other types of waves affect the cascade 
only marginally and the GS95 scaling is expected 
to persist in the compressible medium. Moreover, 
as the Alfven waves are incompressible, their cas- 
cade persists even for supersonic turbulence. In low 
(3 regime the slow modes are sound-type perturba- 
tions moving along magnetic fields with velocity a 
(or, a cos 8 in the direction of wave vector k), where 
a (= y/jPgas/p) is the sound velocity and 9 is the 
angle between the wavevector and magnetic field. In 
magnetically dominated environments a <C Va, the 
gaseous perturbations are essentially static and the 
magnetic field mixing motions are expected to mix 
density perturbations as if they were passive scalar. 
As the passive scalar shows the same scaling as the 
velocity field of the inducing turbulent motions, the 
slow waves are expected to demonstrate GS95 scal- 
ings. The fast waves in low (3 regime propagate at Va 
irrespectively of the magnetic field direction. Thus 
the mixing motions induced by Alfven waves should 
affect the fast wave cascade only marginally. The 
latter cascade should be analogous to the acoustic 
wave cascade and be isotropic. 

To separate the fast, Alfven and slow modes of 
MHD turbulence was the task solved in CL02. Ear- 
lier researchers were probably not much motivated 
to attempt this, as it was generally believed that the 
different modes are messed up in turbulence anyhow. 
To separate the modes, let us consider the displace- 
ment vector £. For example, £/(k) denotes the di- 
rection of displacement for the fast mode with wave 
vector k. In CL02, we showed that fast, slow, and 
Alfven displacement vectors (£/, £ s , and £a, respec- 



MHD TURBULENCE 



3 



1.2 
1 .0 |t 
0.8 
0.6 : 
0.4 : 

0.2 
0.0 





Alfven Modes 



Slow + Fast Modes: 

20 25 



> 



> 

<0 




> 



0.01 




10 15 
t 



I 00 



(B n 2 /4™ r 



<5V 



(B /(4tt Po ) 



1 /2\ 



J / ^"HOi 

Fig. 1. Generation of compressible modes from decaying Alfven Turbulence. Generation is marginal. (Left Panel) 
Time evolution. Initially only Alfven modes exist. As time goes on, compressible modes are spontaneously generated. 

(Middle Panel) Generation of slow modes as a function of Bq. 



Initially, (3 ~ 0.2 and B /^fWpE = 1. 
Generation of fast modes as a function of Bq. From Cho & Lazarian (CL02). 
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Fig. 2. Low [3. Spectra and eddy shapes from driven turbulence with M s ~ 2.2, Ma ~ 0.7, (3 ~ 0.2, and 216 3 
grid points, (a) Spectra of Alfven modes follow a Kolmogorov-like power law. (b) Eddy shapes (contours of same 
second-order structure function, SF2) for velocity of Alfven modes shows anisotropy similar to the GS95 (ry oc r 2 / 3 or 

kn oc k 2 ^ 3 ). The structure functions are measured in directions perpendicular or parallel to the local mean magnetic field 
in real space. We obtain real-space velocity and magnetic fields by inverse Fourier transform of the projected fields, (c) 
Spectra of slow modes also follow a Kolmogorov-like power law. (d) Slow mode velocity shows anisotropy similar to the 
GS95. We obtain contours of equal SF2 directly in real space without going through the projection method, assuming 
slow mode velocity is nearly parallel to local mean magnetic field in low f3 plasmas, (e) Spectra of fast modes are 
compatible with the IK spectrum, (f) The magnetic SF2 of fast modes shows isotropy. We obtain real-space magnetic 
field by inverse Fourier transform of the projected fast magnetic field. Fast mode velocity also shows isotropy. From 
CL02. 

tively) can be written in terms of k|| and k^: 

0( k ) « 1 + a + ^P fc k + k ± k ±7 (1) 
1 + a + V D 



£ s (k) oc fciiki, + 



k ± -k ± , (2) 

-1 + a — VD 

£ A (k) cx kn x kj_, (3) 
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Fig. 3. (a) Iso-contours of equal (velocity) correlation 
from a simulation. The contours represent shape of dif- 
ferent size eddies. The smaller contours (or, eddies) are 
more elongated, (b) Contours generated from analytical 
formula. From Cho, Lazarian, & Vishniac (2002a). 

where k|| (or kj_) is the wave vector parallel (or per- 
pendicular) to the mean field Bo, a = a 2 jV\ = 
07/2, and D = (1 + a) 2 - Aacos 2 9. After 
proper normalization, we get a set of unit bases 
(£/(k),£ s (k),£yi(k)). These unit vectors are well- 
defined and mutually perpendicular except on the 
k|| axis. We can get fast, slow, or Alfven component 
velocity by projecting the Fourier velocity compo- 
nent v(k) onto these unit bases. See CL02 and Cho 
& Lazarian (preprint) for details. 

Results of the CL02 study are shown in Fig p|. 
The calculations confirm the theoretical considera- 
tions provided above. Indeed, both the spectra and 
anisotropics obtained for Alfven and slow modes look 
remarkably similar to the earlier incompressible runs 
(see Fig || where both the data and analytical fit for 
the isocontours of equal correlation are shown) . 

Our more recent results for high (3 plasma are 
shown in Fig |J. They confirm theoretical consider- 
ations in Lithwick & Goldreich (2001). This means 
that for the first time we have a theoretical insight 
and simple scaling relations working for high (3 and 
low f3 MHD turbulence. What will happen in the 
regime of moderate [31 

We would speculate that due to its rapid cas- 
cade the Alfven waves should preserve GS95 scaling. 
Moreover, the passive nature of the slow modes in 
both high and low (3 regimes induce us to think that 
they would mimic the spectrum of Alfven modes in 
the (3 ~ 1 case as well. The decoupling of the fast 
waves and their isotropy in this regime are also likely. 
Thus we may hope that we have a picture of MHD 
turbulence in general and the difference between the 
aforementioned modes would stem from the differ- 
ence in the partition of energy between magnetic and 
gas energies in the compressional modes as f3 varies. 
For instance, most energy of the slow modes in low (3 
plasma is in gas compression, while in high (3 plasma 
the slow modes energy is mostly magnetic. 

How important is the strength of the regular 
magnetic field? In our computations the magnetic 
field was taken to be sufficiently strong. If the mag- 



netic field is weak and strongly entangled by turbu- 
lence, the decoupling of modes is difficult. One may 
wonder to what extend our insight above is appli- 
cable. We believe that if the magnetic field energy 
B 2 /8ir is much smaller than energy pV 2 /2 at a scale 
Z, the magnetic field would not change the dynamics 
of eddies at this scale and therefore the turbulence 
would be hydrodynamic. At smaller scales, however, 
as the energy of eddies decreases (e.g. as I 2 / 3 for 
the Kolmogorov cascade) the magnetic field becomes 
dynamically important with the implication that the 
scaling relations tested for the strong external mag- 
netic field should be applicable. 

3. ION-NEUTRAL DAMPING: A NEW REGIME 
OF TURBULENCE 

In hydrodynamic turbulence viscosity sets a min- 
imal scale for motion, with an exponential suppres- 
sion of motion on smaller scales. Below the viscous 
cutoff the kinetic energy contained in a wavenum- 
ber band is dissipated at that scale, instead of being 
transferred to smaller scales. This means the end of 
the hydrodynamic cascade, but in MHD turbulence 
this is not the end of magnetic structure evolution. 
For viscosity much larger than resistivity, v 3> »/, 
there will be a broad range of scales where viscosity 
is important but resistivity is not. On these scales 
magnetic field structures will be created by the shear 
from non-damped turbulent motions, which amounts 
essentially to the shear from the smallest undamped 
scales. The created magnetic structures would evolve 
through generating small scale motions. As a result, 
we expect a power-law tail in the energy distribu- 
tion, rather than an exponential cutoff. This com- 
pletely new regime for MHD turbulence was reported 
in Cho, Lazarian & Vishniac (2002c). Further re- 
search showed that there is a smooth connection be- 
tween this regime and small scale turbulent dynamo 
in high Prandtl number fluids (see Schekochihin et 
al. 2002). 

In partially ionized gas neutrals produce viscous 
damping of turbulent motions. In the Cold Neutral 
Medium (see Draine & Lazarian 1999 for a list of the 
idealized phases) this produces damping on the scale 
of a fraction of a parsec. The magnetic diffusion in 
those circumstances is still negligible and exerts an 
influence only at the much smaller scales, ~ lOOfcm. 
Therefore, there is a large range of scales where the 
physics of the turbulent cascade is very different from 
the GS95 picture. 

Cho, Lazarian, & Vishniac (2002c) explored this 
regime numerically with a grid of 384 3 and a physical 
viscosity for velocity damping. The kinetic Reynolds 



MHD TURBULENCE 



■5 




Fig. 4. High (3. Spectra and eddy shapes from driven turbulence with M s ~ 0.5, Ma ~ 0.7, (3 ~ 4, and 216 3 grid 
points. See captions in Fig. 0. From Cho & Lazarian (in prep.). 
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Fig. 5. New Regime of MHD turbulence: viscously 
damped, (a) (Left) Incompressible case. A new iner- 
tial range emerges below the viscous cut-off at k ~ 7. 
From Cho, Lazarian, & Vishniac (2002c). (b) (Right) 
Due to numerical reasons the inertial range is truncated 
when we use our compressible code. However, the new 
regime persists and causes density fluctuations. From an 
upcomming paper. 

number was around 100. The result is presented in 
Fig. |a. 

A theoretical model for this new regime and its 
consequences for stochastic reconnection (Lazarian 
& Vishniac 1999) can be found in Lazarian, Vish- 
niac, & Cho (2002). It explains the spectrum fc _1 
as a cascade of magnetic energy to small scales un- 
der the influence of shear at the marginally damped 
scales. Moreover, our work suggests that the mag- 



netic fluctuations protrude to the decoupling scales 
and cause the renewal of the MHD cascade there. 
Earlier work, e.g. Lithwick & Goldreich (2001) ar- 
gued that the turbulent cascade survives ion-neutral 
damping only when a high degree of ionization is 
present. In view of the finding a revision of a few 
earlier theoretical conclusions is necessary. It should 
be noted that our conclusion about the resumption 
of turbulence at small scales is consistent with obser- 
vations (Spanglcr 1991; 1999) that do not show any 
change of the observed electron scintillation spec- 
trum at the ambipolar damping scale. 

We show our results for the compressible fluid 
in Fig ||b. The inertial range is much smaller due 
to numerical reasons, but it is clear that the new 
regime of MHD turbulence persists. The magnetic 
fluctuations, however, compress the gas and thus 
cause fluctuations in density. This is a new (al- 
though expected) phenomenon compared to our ear- 
lier incompressible calculations. These density fluc- 
tuations may have important consequences for the 
small scale structure of the ISM. We may specu- 
late that they might have some relation to the tiny- 
scale atomic structures (TSAS). Heiles (1997) intro- 
duced the term TSAS for the mysterious H I absorb- 
ing structures on scales from thousands to tens of 
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AU, discovered by Dieter, Welch & Romney (1976). 
Analogs are observed in Nal and Call (Meyer & 
Blades 1996; Faison & Goss 2001; Andrews, Meyer, 
& Lauroesch 2001) and in molecular gas (Marscher, 
Moore, & Bania 1993). 

Our calculations are applicable on scales from the 
viscous damping scale (determined by equating the 
energy transfer rate with the viscous damping rate; 
~ 0.1 pc in the Warm Neutral Medium with n = 

0. 4 cm' 3 , T= 6000 K) to the ion-neutral decoupling 
scale (the scale at which viscous drag on ions be- 
comes comparable to the neutral drag; « 0.1 pc). 
Below the viscous scale the fluctuations of magnetic 
field obey the damped regime shown in Figure |5|b 
and produce density fluctuations. For typical Cold 
Neutral Medium gas, the scale of neutral-ion decou- 
pling decreases to ~ 70 AU, and is less for denser 
gas. TSAS may be created by strongly nonlinear 
MHD turbulence! 

4. DISCUSSION 

In this paper we have discussed the new outlook 
onto compressible MHD turbulence. Contrary to 
common beliefs the compressible MHD turbulence 
does not present a complete mess, but demonstrates 
nice scaling relations for its modes. A peculiar fea- 
ture is that those relations should be studied locally, 

1. e. in the frame related to the local magnetic field. 
However, such a system of reference is natural for 
many phenomena, e.g. for cosmic ray propagation. 
Recent application of the scalings obtained for com- 
pressible turbulence have shown that fundamental 
revisions are necessary for the field of high energy 
astrophysics. For instance, Yan & Lazarian (2002) 
demonstrated that fast modes dominate cosmic ray 
scattering even in spite of the fact that they are sub- 
jected to collisional and collisionless damping. This 
entails consequences for models of cosmic ray prop- 
agation, acceleration, elemental abundances etc. 

Advances in understanding of MHD turbulence 
have very broad astrophysical implications. The 
fields affected span from accretion disks and stars 
to the ISM and the intergalactic medium in clusters. 
Turbulence is known to hold the key to many astro- 
physical processes. It was considered too messy by 
many researchers who consciously or subconsciously 
tried to avoid dealing with it. Others, more brave 
types, used Kolmogorov scalings for compressible 
strongly magnetized gas, although they did under- 
stand that those relations could not be true. 



Recent research in the field provides the scaling rela- 
tions and insights that will contribute to many areas 
of research. 

AL and JC acknowledge support by the NSF 
grant AST0125544. 
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